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Abstract 

We argue that study of the cross section of coherent photo(electro) production of vector mesons 
off nuclear targets provides an effective method to probe the leading twist hard QCD regimes of color 
transparency and perturbative color opacity as well as the onset of black body limit (BBL) in the soft and 
hard QCD interactions. In the case of intermediate energies we use the Generalized Vector Dominance 
Model (GVDM) to take into account coherence effects for two distinctive limits - the soft interactions for 
production of p and p'-mesons and the color transparency regime for production of charmonium states. 
We demonstrate that GVDM describes very well p-meson coherent photoproduction at 6 < E~f < 10 
GeV and predict an oscillating energy dependence for the coherent charmonium production. In the limit 
of small X we find that hard QCD leads to onset of the perturbative color opacity even for production 
of very small onium states, like T. The advantages of the process of coherent dijet photoproduction and 
hard diffractive processes in general for probing the onset of BBL and measuring the light-cone wave 
function of the photon in a hard scattering regime where decomposition over twists becomes inapplicable 
are explained. We apply this analysis to the study of the photon induced coherent processes in Ultra 
Peripheral Collisions of ions at LHC and demonstrate that the counting rates will be sufficient to study 
the physics of color opacity and color transparency at the energies beyond the reach of the electron- 
nucleon(nucleus) colliders. 

1 INTRODUCTION 



It appears that in the next decade the photon - nucleus interactions will be in the forefront of the small 
X QCD dynamics studies. This is due to the possibilities of the studying coherent (and for some channels 
incoherent) photon - nucleus interactions at energies which exceed at least by a factor of 10 the energies of 
electron-nucleon interactions at HERA. Such studies will be feasible at LHC within the program of the study 
of the Ultra Peripheral Collisions (UPC) m^lEl- This opens a challenging opportunity to get answers to a 
number of the fundamental questions which could be investigated in the coherent processes: how to achieve 
new QCD regime of strong interaction with small coupling constant, how interactions depend on the type of 
the projectile and how they change with an increase of the size/thickness of the target, etc. Several regimes 
appear possible depending on the incident energy and the target thickness. A hadronic projectile (proton, 
pion, etc) high-energy interactions with the nucleus rather rapidly approach the black body limit (BBL) 
in which the total cross section of the interaction is equal to 2'n:R\ {Ra — 1.2 • A^^^). Another extreme 
limit is the interaction of small size projectiles (or wave packages). In this case in a wide range of high 
energies the system remains almost frozen during the passage through the nucleus and the regime of color 
transparency is reached in which the interaction of the small size projectile with a nucleus is rather weak and 
proportional to A. This color transparency phenomenon in hard high energy projectile-nucleus interactions 
has been recently observed experimentally in the exclusive production of two jets by pions in the coherent 
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diffraction off nuclei. At higher energies interactions of sniaU dipole is expected to reach the regime of the 
perturbative color opacity. In this regime the small size projectile still couples to the gluon field of the target 
via the (skewed) gluon density of the target, Ga{x, Q^), like in the color transparency regime. However, the 
scattering amplitude is not cx A due to the leading twist (LT) shadowing, Ga/AGn < 1- The onset of gluon 
shadowing tames somewhat the increase of Ga{x, Q^) between x ^ 10~^ and x ^ 10^^ and slows down the 
increase of the small dipole - nucleus interaction with energy. The resulting taming is not strong enough to 
prevent the leading twist approximation result for the total inelastic cross section from reaching and even 
exceeding the BBL. This violation of unitarity is an unambiguous signal that for sufficiently small x the LT 
approximation breaks down. 

An important practical issue for the studies of the small x dynamics is whether the interactions of small 
dipoles with several nucleons of the target are strongly modified by the leading twist shadowing dynamics 
0] or one can neglect the LT gluon shadowing, like in the model of McLerran and Venugopalan^, and focus 
on the higher twist effects which are often modeled in the impact parameter space eikonal model 6 . If 
the leading twist shadowing was small and only higher twist effects were taming the increase of the dipole 
- nucleus cross section, the break down of the DGLAP approximation would occur at rather large x. The 
break down of DGLAP may result in onset of the BBL or taming of the cross sections at smaller values. We 
will argue below that the relative importance of leading and higher twist interactions could be experimentally 
resolved using coherent onium photoproduction. 

Note also that if the LT gluon shadowing effects were small enough, the BBL for the interaction of qq 
dipoles of the size > 0.3 — 0.4/r7i could be reached for central collisions with heavy nuclei already at a; > 10^'^ 
that is in the kinematics where In x effects in the evolution of the parton densities are small. In any case, 
whatever is the limiting behavior for the interaction of the small size dipoles with heavy nuclei, it is of major 
theoretical interest since it represents a new regime of interactions when the leading twist approximation 
and therefore the whole notion of the parton distributions becomes inapplicable for the description of hard 
QCD processes in the small x regime. It is worth emphasizing that on the top of providing higher parton 
density targets, nuclei have another important advantage as compared to the nucleon target. It is a weak 
dependence of the scattering amplitude on the impact parameter for a wide range of the impact parameters 
(in fact, one can combine light and heavy nuclei to study the dependence of the amplitudes on the nuclear 
thickness) . On the contrary, in the nucleon case the scattering at large impact parameters is important even 
at very small x working to mask the change of the regime of the interactions at small impact parameters. 

Theoretical studies of the limiting behavior of the small dipole - heavy nucleus cross sections did not lead 
so far to definitive results. It is conceivable that the QCD dynamics will stop the increase of the cross section 
at central impact parameters at the values significantly smaller than allowed by the BBL. In the following 
discussion to emphasize the qualitative difference of the new regime we will use for simplicity the extreme 
hypothesis that the impact factor at small impact parameters reaches the value of one corresponding to the 
BBL pattern of the elastic and inelastic cross sections being equal. 

There exist many processes where the projectile wave function is a superposition of configurations of 
different sizes, leading to the fluctuations of the interaction strength. In this respect, interactions of real 
and virtual photons with heavy nuclei provide unique opportunities since the photon wave function contains 
both the hadron-like configurations (vector meson dominance) and the direct photon configurations (small 
qq components, heavy quark-antiquark components). The important advantage of the photon is that at high 
energies the BBL is manifested in diffraction into a multitude of the hadronic final states (elastic diffraction 
7 — > 7 is negligible) while in the hadron case only elastic diffraction survives in the BBL, and details of the 
dynamics responsible for this regime remain hidden. Moreover, one can post-select a small size or a large 
size initial state by selecting a particular final state, containing for example cc 66, or a leading light meson. 
Such post-selection is much easier in the photon case than in the hadron case since the distribution over the 
sizes of the configurations is much broader in the photon case. 

Spectacular manifestations of BBL in (virtual) photon diffraction include strong enhancement of the 
large mass tail of the diffractive spectrum as compared to the expectations of the triple Pomeron limit, large 
cross section of the highpt dijet production 0. We emphasize that the study of the diffractive channels will 
allow to distinguish between two extreme scenarios: large suppression of the cross section due to the leading 
twist shadowing and nonlinear regime of the BBL. Investigation of the coherent diffraction in BBL would 
allow to perform unique measurements of various components of the light cone wave function of the photon, 
providing a much more detailed information than similar measurements in the regime where leading twist 
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dominates. 

In this review we summarize our recent studies of the various regimes of the coherent photoproduction 
off nuclei |7ll5ll^ [TUl IllL IT^ I13L I14| : the onset of the BBL regime, phenomenon of color transparency and 
perturbative color opacity related to the leading twist nuclear gluon shadowing, and the pattern of soft QCD 
phenomena in the proximity to the black body limit. In addition we consider hard leading twist diffraction 
off nuclei in DIS. We also outline how these effects can be studied in UPC collisions and provide a comparison 
with the first UPC data from RHIC. 



2 VECTOR MESON PRODUCTION OFF NUCLEI IN THE 
GENERALIZED VECTOR DOMINANCE MODEL AT IN- 
TERMEDIATE ENERGIES 

2.1 Outline of the model 

Our consideration of the coherent photoproduction processes of light flavor and hidden charm is based on the 
use of the eikonal approximation that is the Glauber model modified to take into account finite longitudinal 
momentum transfer The photoproduction cross section a-yA^pAi^^) is given in the Glauber model by 
the general expression 
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(1) 



Here lo^ is the photon energy, /cyis the vector meson momentum, = t± — — t, —tmin — t-t is the 

longitudinal momentum transfer in the j — V transition, and T{b) is the diffractive nuclear profile function. 
Depending on the considered vector meson production process, the Glauber approach can be combined with 
either the traditional vector dominance model (VDM) (for the detailed review see Ref. P^) or with the 
generalized vector dominance model (GVDM) |17l 1181 IT?^ . More properly this latter approximation should 
be called the Gribov-Glauber model [20] because the space-time evolution of high energy processes is different 
in quantum mechanical models and in quantum field theory, and therefore theoretical foundations of the 
high-energy model are different. In particular, in QCD in difference from quantum mechanics a high-energy 
projectile interacts with all nucleons at the same impact parameter almost at the same time [23 . Due to the 
cancellations between diagonal and nondiagonal transitions the GVDM allows to take into account the QCD 
effect of the suppression of interaction of spatially small quark-gluon wave packages with a hadron target - 
the so called color transparency phenomenon. Namely, cancellations occur between the amplitudes of the 
photon transition into a vector state Vi with subsequent conversion to a state V2 and direct production of the 
state V2, etc. The importance of the nondiagonal transitions reveals itself in the precocious Bjorken scaling 
for moderately small x ^ 10^^ as due to the presence in the virtual photon of hadron-like and point-like type 
configurations |21) . These amplitudes are also crucial for ensuring a quantitative matching with perturbative 
QCD regime for < few GeV'^ ,22,. The amplitude of the vector meson production off a nucleon can be 
written within the GVDM as 

v4(7 + TV ^ Vj +N) = Y -^A{Y, + ^ + TV), (2) 

where fvi are expressed through r{Vi — > e+e"). Calculation of the vector meson production amplitude off 
nuclei within the Glauber approximation requires taking into account both the nondiagonal transitions due 
to the transition of the photon to a different meson V' in the vertex 7 V^' and due to a change of the 
meson in multiple rescatterings like j —^ V —f V' V . This physics is equivalent to inelastic shadowing 
phenomenon familiar from hadron- nucleus scattering [20]. Then in the optical limit {A^ 1) of the Glauber 
multistep production theory one can introduce the eikonal functions , ^v.v {b, z), which describe propagation 
of the produced objects through the medium and are related to the diffractive profile function by expression 

r(6) = lim (3) 
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Within the optical limit of the Glauber based GVDM with accuracy 0{y/a^^) the eikonal functions ^v,V' {b, z) 
are determined by the solutions of the coupled channel equations 



dz 



2ikv 



UjA^VAib,z)e' 



E UvA^v,A{b, z)e'^^^""'<^v' (6, z) 



(4) 



with the initial condition <i>v'.v''(6, — cx)) — 0. The exponential factors ex'p\iq^^ z\ are responsible for the 

coherent length effect, «,j = ^,V,V', q\^-' = — '-^^^ — -. The generalized Glauber -based optical potentials in 
the short-range approximation are given by the expression 

UiA^jA{b,z) = -47r/ijv^jAr(0)£'(fo, z). 



(5) 



Here fiN^jN{0) are the forward elementary amplitudes, and g(6, z) is the nuclear density normalized by the 
condition J dbdz g{b, z) = A. We calculated g{b, z) in the Hartree-Fock-Skyrme (HFS) model which provided 
a very good(with an accuracy « 2%) description of the global nuclear properties of spherical nuclei along 
the periodical table from carbon to uranium|23j and the shell momentum distributions in the high energy 
(p,2p)[2^ and (e,e'p)j^ reactions. If the nondiagonal rescattering amplitudes fvN^VN = 0, one can easily 
integrate Eq^ and using the expression for the elementary amplitude 



fvN-*VN{t — 0) 
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47r 



obtain the expression for the photoproduction cross section: 
da-^iA^VA _ da^N^VN{t = 0) 



dt 



dt 



d^bdze'^'-''pib,z)e""'' -e 



1 _t 
' 2 " ViV 



^(s)(l-mv7v) J p(b,z')dz' 



(6) 



well known from early seventies (see for example [Ifcip. 



2.2 Production of light vector mesons 

We have used the GVDM to describe coherent photoproduction of hadronic states of M < 2 GeV off nuclei 
and consider the onset of BBL in the soft regime^ . In Ref . j26j the simplest nondiagonal model was considered 
with two states p and p' . Then the GVDM comprises elementary amplitudes 

f-jN^pN — -rfpN^pN + -^fp'N^pN, 
JP Jp' 

6 6 

flN^p'N = -7—fp'N^p'N + -rfpN^p'N- (7) 
Jp' Jp 

It was assumed that both p and p' have the same diagonal amplitudes of scattering off a nucleon. The ratio 
of coupling constants was fixed: fpi / fp = \/3, while the ratio of the nondiagonal and diagonal amplitudes 
A{p + N ^ p' + N)/A{p + N ^ p + N) = — e, and the value cr^'jj were found from the fit to the forward 
J + A ^ p + A cross sections measured at oj^ =6.1, 6.6 and 8.8 GeV 27 . One should emphasize here that 
in such GVDM extension p'-meson approximates the hadron production in the interval of hadron masses 
AM^ ^ 2 GeV^. Thus the values of the production cross section refer to the corresponding mass interval. 

We refined this model in The dependence on the nuclear structure parameters was diminished by 
calculating the nuclear densities in the Hartree-Fock-Skyrme (HFS) approach. Next, we used in all our 
calculations the parameterization of P2] for the pN amplitude which was obtained from the fit to the 
experimental data on photoproduction off the proton target. The value of e was fixed at 0.18 to ensure 
the best fit of the measured differential cross section of the p-meson photoproduction off lead at oj^ — 6.2 

^In our calculation we neglect the triple Pomcron contribution which is present at high energies. This contribution though 
noticeable for the scattering off the lightest nuclei becomes a very small correction for the scattering of heavy nuclei due the 
strongly absorptive nature of the interaction at the central impact parameters. 
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GeV and t± — 0.001 GeV^. Note that this value of e leads to a suppression of the differential cross 
section of the p-photoproduction in7+p^ p + pbya factor of (1 — w 0.80 practically coinciding 

with phenomenological renormalization factor R = 0.84 introduced in 28 to achieve the best fit of the 
elementary p-meson photoproduction forward cross section in the VDM which neglects mixing effects. With 
all parameters fixed we calculated the differential cross sections of p-production off nuclei and found a good 
agreement with the data p7, see a detailed comparison in 9|. In view of a good agreement of the model 
with the data on p-meson production in the low energy domain we used this model to consider the p-meson 
photoproduction at higher energies of photons. The increase of the coherence length with the photon energy 
leads to a qualitative difference in the energy dependence of the coherent vector meson production off light 
and heavy nuclei and to a change of the A-dependence for the ratio of the forward p' and p-meson production 
cross sections between uj^ ~ 10 GeV and u-^ ^ 40 GeV (Fignj. The observed pattern reflects the difference 
of the coherence lengths of the p-meson and a heavier p'-meson which is important for the intermediate 
photon energies < 30 GeV. The corrections due to nondiagonal transitions are relatively small 15%) for 
the case of p production off a nuclei. As a result, we find that the GVDM cross section is close to the one 
calculated in the VDM for heavy nuclei as well. Situation for p' production is much more interesting. The 
cross section of p' production off a nucleon is strongly suppressed as compared to the case when the p <-> p' 
transitions are switched off. The extra suppression factor is w 0.5. The non-diagonal transitions disappear 
in the limit of large A (corresponding to the BBL) due to the condition of orthogonality of hadronic wave 
functions in accordance with the general argument of Gribov |17| . Therefore in the limit oi A —^ oo we 
expect that the relation 



da{"/ + hi+ A)/dt _ cr(e+e" ^ /ii) 

da{-f + A^h2 + A)/dt |^^^ ~ cr(e+e- /12) 



{12/ hf, (8) 



should be fulfilled for the productions of states /ii, /12 of invariant masses , A/| at t± ~ 0. Indeed we have 
found from calculations that in the case of the coherent photoproduction off lead the nondiagonal transitions 
becomes strongly suppressed with increase of the photon energy. As a result the p'/p ratio increases, exceeds 
the ratio of the 7p Vp forward cross sections calculated with accounting for p ^ p' transitions already 
at uj^ > 50 GeV and becomes close to the value of fp/fpi which can be considered as the limit when one 
can treat the interaction with the heavy nucleus as a black one. It is worth noting here that presence 
of nondiagonal transitions in terms of the formalism of the scattering eigen states [21] corresponds to the 
fluctuations of the values of the interaction cross sections for the real photon. The GVDM discussed in the 
paper leads to small (ex 10%) color transparency effects at intermediate energies for the cross section of 
semiinclusive photoproduction processes. Really, this model corresponds to the propagation of states with 
cross sections: w a{VN){l ± e) . In the case of electroproduction e should be significantly larger: 



e^lfL= / r(p^e+e-) ^, r(p'^e+e-) 
fp y mp mp> 

This equation follows from Eq. Q where the left hand side is put to zero because the cross section of the elastic 
vector meson electroproduction rapidly decreases with Q^. The presence of the CT phenomenon within the 
GVDM leads to a substantial modification of the pattern of the approach to BBL. The nondiagonal transitions 
become more important with the increase of Q^, leading to an enhancement of the effects discussed above. 
In particular, the fluctuations of strengths of interaction would lead at large to the color transparency 
phenomenon. Presence in GVDM of significantly different masses of p, p', ... makes it impossible to describe 
all fluctuations of strengths of interaction in terms of one coherent length. 



2.3 Photoproduction of J /ip,ip' -mesons at intermediate energies. 

There is a qualitative difference between GVDM description of the light meson production described above 
and production of charmonium states. In the case of the photoproduction of light mesons the soft physics 
dominates both at low and high energies. Hence, the Gribov-Glauber model should be applicable in a wide 
range of energies. At the same time the space-time picture of the process changes with energy. At low 
energies the meson is formed at the distances smaller than the typical interaction length of a meson in the 
nucleus. With the increase of energy the formation length starts to exceed the nuclear size and a photon 



5 




Figure 1: The energy dependence of the p-photoproduction cross section and the A-dependence of the p' j p 
photoproduction cross sections calculated in the GVDM+Glauber model. 
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converts to a system of qq and higher Fock-components before the target, and one has to account for the 
interactions of this system with the media and subsequent transition of the system to a vector meson. 

The situation is much more involved in the case of (charm)onium production. At high energies both 
the coherence length Ic « 2u>-ymy^ and the formation length If « 2a;^[m^, — Wj^^] (distance on which 
the squeezed qq pair transforms into the ordinary meson) are large and the color transparency phenomenon 
reveals itself. In particular, it explains the fast increase of the cross section with energy observed at HERA (for 
reviews of small x phenomena see Ref. [HOIE]). Consequently, the value of the cross section extracted from 
the charmonium photoproduction characterizes the interaction of the squeezed cc pair with a nucleon rather 
than the charmonium-nucleon interaction. In the high energy limit the very small interquark distances in the 
wave function of the photon dominate, and one has to treat the interaction of small dipoles with nuclei. In 
this case the eikonal approximation gives a qualitatively wrong answer since it does not take into account the 
leading twist effect of the gluon shadowing(see the detailed analysis in while the leading twist analysis 
predicts large shadowing effects [S], see discussion in section 3. On the other hand, at the intermediate 
energies when onium states are formed inside the nucleus the nonperturbative effects at a transverse distance 
scale comparable to the charmonium size becomes important. The hadronic basis description would be more 
relevant in this case. However the VDM which takes into account only diagonal vector meson transitions 
does not account properly for the basic QCD dynamics of interaction. In particular, SLAC data j^Hl show 
that (7-y+N~>4>'N ~ 0.15 • a^+N^ipN which within VDM corresponds to: ly-ip' n / <y j /ipN ~ 0-7. This conclusion 
is in evident contradiction with the QCD expectation that the hadron interaction cross section should be 
scaled approximately as the transverse area occupied by color: cr^'jv/fj/)/,iv oc Thus in this naive 

QCD picture cr{4>'N)/a{J/4>N) ~ 4. A QCD explanation of such failure of VDM is based on the observation 
|34[I35| that in photoproduction of both the J/'f/; and tp' mesons the small relative distances ^ l/rric dominate 
in the cc component of the photon wave function. Therefore the cross sections which enter into the ratio 
of J/ip and ip' yields are cross sections of the interaction of the small dipoles, not genuine mesons. The 
suppression of the production of ijj' in this picture is primarily due to a smaller leptonic decay width of tp' (a 
factor of 1/3). A significant additional suppression comes from importance of more massive cc intermediate 
states in the photon wave function (> Mj/^^M^i respectively). The exact value of the suppression is 
sensitive to the details of the onium wave functions and the dependence of the dipole - nucleon cross section 
on the size of the dipole 36, 37 . 

The dominance of small cc configurations in photoproduction processes is relevant for the significant 
probability of nondiagonal J /ip «-> -0' diffractive transitions. The GVDM adjusted to account for the color 
screening phenomenon 34, 35 allows to take into account QCD dynamics using a hadronic basis. Such 
a description is limited to the regime of small coherence lengths (medium energies), where leading twist 
shadowing is not important. We use the GVDM to consider the coherent photoproduction of hidden charm 
mesons off nuclei at moderate photon energies 20 GeV < w-^ < 60 GeV where the coherence length for 
the transition 1^ is still close enough to the internucleon distance in nuclei while the formation length 
If is comparable to the radii of heavy nuclei. In this energy range produced charmonium states have 
a noticeable probability to rescatter in sufficiently heavy nuclei. Therefore, one would be able to reveal 
the fluctuation of the charmonium-nucleon interaction strength as due to the diagonal -ipN ipN and 
nondiagonal '^pN <S4> ip' N rescatterings for moderate energies. The GVDM which we outlined above takes 
into account the coherence length effects via the Glauber model approximation. The key distinction from the 
/5— meson case is choosing the parameters of the model to account for the space-time evolution of spatially 
small cc pair. It is important that the inelastic shadowing corrections related to the production of higher 
mass states |2()j are still insignificant. A reasonable starting approximation to evaluate the amplitude of the 
charmonium-nucleon interaction is to restrict ourselves to the basis of J/ ip and tp' states for the photon wave 
function in EqEl Then, similarly to the considered above case of p, p' production we have two equations 

comprising six elementary amplitudes. The charmonium-nucleon coupling constants: = 10.5 ± 0.7, 

/^ 

and = 30.9 ± 2.6 are determined from the widths of the vector meson decays V — s- ee. Since in 
the photoproduction processes cc pair is produced within the spatially small configuration one can neglect 
for a moment by the direct photoproduction amplitude and obtain the approximative relations between 
rescattering amplitudes 
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fj/ipN- 



■J/ipN- 



(9) 



Large values for nondiagonal amplitudes f^i n~*j/^n ~ '^^■'^fj/ipN—>j/^N £^re a characteristic QCD property 
of hidden charm and beauty meson-nucleon interaction. Note that the negative sign of the nondiagonal 
amplitude is dictated by the QCD factorization theorem. A positive sign of the forward photoproduction 
f.yN^j/4}{')p')N amplitudes as well as the signs of the coupling constants fj/^ and are determined by the 
signs of the charmonium wave functions at r=0. From the approximative estimates above it also follows 
that (T^'N ~ 9(Tj/^7v This is much larger than (t^'n ^ 20 mb suggested by analyses of the data on ijj' 
absorption in nucleus- nucleus collisions ,38. .39) . When combined with the SLAC data 001) this corresponds 
to (JtP'n I'^j/xpN ~ 5 6 with large experimental and theoretical errors. 

To fix elementary amplitudes more accurately within GVDM (see ^2 fo^' details) we parameterized the 
elementary photoproduction cross section in the form used by the experimentalists of HERA to describe 
their data. This form has no firm theoretical justification but it is convenient for the fit 



0.4 



(10) 



Here s — 2LUjm]\j is the invariant energy for the photon scattering off a free nuclcon, sq = 40 GeV^ is the 
reference point and i^2g (t) is the two-gluon form factor of a nucleon. 

We also used the SLAC results for the forward 77V — > J /-ipN cross section and for the ratio 



dfj. 



dt 



0.15 



dt 



(11) 



Besides, we fixed the J/tpN cross section aj/^jy ~ (3.5 ± 0.8) mb as measured at SLAC^OI and used the 
reasonable assumption about the energy dependence of this cross section as the sum of soft and hard physics: 



3.2 mb I — 



0.08 



0.3 mb I — 

■So 



0.2 



(12) 



The existence of a hard part of the J/ipN cross section is consistent with the GVDM, because the pho- 
toproduction amplitude has a stronger energy dependence than the Pomeron exchange, i.e. soft scattering 
amplitudes. Next, we used this input, the optical theorem and the well known Gribov-Migdal relation 



STT d 3/, 



2 d\ns 



(13) 



to determine from Eq|2] all elementary amplitudes in the discussed energy range. In particular, we found 
the value ct^'at ~ 8 mb. However, the input parameters of the GVDM, namely the experimental cross 
sections of the forward elementary photoproduction and, especially, the value of (Jj/tpNi are known with 
large uncertainties. In result, we obtained elementary cross sections changing in the ranges: 2.5 mb < 
cTj/^N < 5 rnfo and 6mb < ct^'jv < 12 mb, in the discussed energy range. We checked how a variation of 
aj/^N within the experimental errors influences our results. 

In our papers^^liCSl we analyzed the coherent photoproduction of charmonia off light(Si) and heavy(Pb) 
nuclei. Here we present the cross sections of the coherent photoproduction of J/tp and ip' off Ca. The energy 
dependence of these cross sections is compared (Fig|2Ji) to that obtained in the Impulse Approximation 
where all rescatterings of the produced vector mesons are neglected, and the cross section is given by the 
simple formula 



d(T^A^VA{s,t) _ da^N^VN{s,tmin) 



dt 



dt 



"db 



g{b, z) 



(14) 



The distinctive feature of the coherent charmonium photoproduction is that differential cross sections 
oscillate as a function of the photon energy. The major source for such a behavior at intermediate energies 
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Figure 2: The coherent photoproduction of charmonia off Ca: a). Energy dependence of photoproduction 
cross sections; b). Ratio of cross sections calculated within the GVDM and lA for J/tjj', c). The same for 
d). Ration of the forward ■(/''-production cross section to that for J/tp. The filled areas show the variation of 
the cross sections due to the uncertainty of the experimental J/i/iN cross section. 
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Figure 3: The A-dependence of the coherent charmonium photoproduction cross sections calculated in the 
Generalized Vector Dominance Model for different energies of the photon. 
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is oscillating behavior of the longitudinal nuclear form factor at the relatively large values of imm in the 
photoproduction vertex. 

The cross section of the J j 4> photoproduction in the GVDM is close to that calculated in the Impulse 
Approximation - the shapes of curves are very similar and the values of the cross sections are only slightly 
reduced at energies below 40 GeV(Fig|2b)- This indicates that it would be very difficult to extract the total 
J/tpN cross section from such a measurement. At the same time, since the the nuclear form factor is known 
with a high precision from the high energy elastic electron-nucleus scattering, one can use the coherent 
J/tp photoproduction off the spherical nuclei to determine the elementary J/ipN amplitude in a wide 

range of energies for t ~ 0. 

The picture is qualitatively different for the coherent ip' photoproduction off nuclei. Due to the higher 
threshold (i?^ — E'^l^) « 3 GeV, the direct V'' production off nuclei is suppressed in the impulse approxi- 
mation (for the same incident energy) as compared to the J/?/' production by the nuclear form factor. Since 

qj^ w qj'^^^, the minima of the impulse approximation distribution are shifted. The contribution of 

the nondiagonal term 7 J /ip — > ^' essentially increases the ip' yield and produces an additional shift of 
the minima in the spectrum to lower photon energies (Fig. This results in significant effects, especially, 
if one analyses the data as a function of the mass number A at different energies to extract A-dependence 
of the J/ipN cross sections. It is seen from FigElthat the oscillating form factor significantly influences 
the A-dependence that complicates extracting of the genuine J/ipN and ijj' N cross section from such anal- 
ysis. The oscillation of the cross section with energy can be used to measure in a new way the elementary 
charmonium photoproduction amplitudes as well as the charmonium-nucleon amplitudes using the coherent 
charmonium photoproduction off light nuclei [T^. In particular, at the photon energies « O.liRArn^, 
(Fig. the main contribution to the cross section originates from the nondiagonal rescattering. As a 
result, one can extract from the data the nondiagonal elementary J/ipN — > ijj'N amplitude by measuring 
the ratio of the tp' and J/tp yields at zero production angle, because other inputs to this ratio such as the 
elementary •jN J/ipN amplitude are fixed from the J/ip production. Moreover, since other parameters 
enter both in the numerator and denominator, the major uncertainties are canceled out. On the other 
side, the energy dependence of this ratio should originate primarily due to the contribution of the direct ip' 
production. Thus, one would be able to determine the 7iV tp'N amplitude from the measurement of the 
tp'-to-J/tp ratio. We want to emphasize that the suggested procedure for extracting the nondiagonal am- 
plitude and amplitudes of direct J/tp and ip' photoproduction from the nuclear measurements is practically 
model independent. 

It is much more difficult to determine the diagonal J/tpN J/tpN and ip'N ip'N amplitudes, which 
are relevant for the suppression of the charmonium yield in heavy ion collisions. 

A naive diagonal VDM with (Jtot{J /ip^) based on the SLAG data ^Qj leads to a rather significant 
suppression of the J/ipi yield: « 10%-=- 15% for light nuclei, and « 30%-;- 40% for heavy nuclei. However, we 
find jl2j a strong compensation of the suppression due to the contribution of the nondiagonal transitions. 
As a result we find that overall the suppression does not exceed 5% 10% for all nuclei along the periodical 
table. Hence, an extraction of the diagonal amplitudes from the measured cross sections would require 
a comparison of high precision data with very accurate theoretical calculations including the nondiagonal 
transitions. Another strategy is possible if the elementary — > J/ipN, and 7iV ip'N photoproduction 
amplitudes, as well as the nondiagonal amplitude J/ipN <-> ip' N would be reliably determined from the 
medium energy data. In this case it would be possible to determine the imaginary parts of the forward 
diagonal amplitudes from the GVDM equations 10) 



The suggested procedure allows to determine all elementary amplitudes with a reasonable precision from 
the measurements of J/ip and ip' photoproduction off the light nucleus at the medium photon energies. The 
measurement of the coherent charmonium photoproduction off nuclei in this energy region are planned at 
SLAG |41| . The measurements of coherent production off heavy nuclei and of the quasielastic production 
with parameters of GVDM fixed from the analysis of light nuclei would provide a critical test of the model. 




(15) 



(16) 
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The main limitation of the suggested procedure is the restriction of the hadronic basis to the two lowest 
IS, 2S charmonium states with the photon quantum numbers: J /tp , tp' . This is one of the key approximations 
in the discussed approach. It seems quite reasonable since in the coherent production of vector meson states 
off nuclei at moderate energies contribution of high mass states is more suppressed by the target form factor. 
However, we neglected by ip" which is nearly degenerate in mass with -0': AM = 91 MeV. The properties of 
these two states are described well in the charmonium model | 42[I43 |. In this model '(/'"-meson is described as 
a ^Di state with a small admixture of the S wave, while i/j'-meson has a small D-wave admixture. Namely 

\'iP') = cos6'|2S') +sin6l|lD) , 

IV'") = cos6'|li:>) -sin6l|2S') . (17) 

Since only the S-wave contributes to the decay of ip states into e^e^ (at least in the non-relativistic char- 
monium models) the value = 19° ± 2° can be determined from the data on the e+e^ decay widths 
T{tj)' e+e") = 2.19 ± 0.15 KeV and the r(i/'" e+e") = 0.26 ± 0.04 KeV. Due to the smah difference of 
masses between the ■0' and ip" mesons the produced S-wave cc-state does not loose coherence while going 
through the media at any conceivable energies. The soft interactions cannot transform the S'-state to D-state 
with any significant probability. In the soft QCD processes data show that the cross sections of exclusive 
nondiagonal transitions are negligible for the forward angle scattering. The same conclusion is valid in the 
PQCD model for the charm dipole-nucleon interactions. Thus, it is more appropriate to use the \S - 2S 
basis for description of the propagation of cc through the nucleus. The only resulting change is an increase 
of /2s by a factor of 1/ cos 9 as compared to which is within the uncertainties of the model. Since the 
25' state ultimately transforms into tp' and V'" we predict ^]: ^ 

^ = t.A0)^0.1 . (18) 

Influence of the higher mass resonances is expected to be even weaker - the constants 1/fv relevant 
for the transition of a photon to a charmonium state V rapidly decrease with the resonance mass. This 
is because the radius of a bound state, ry, is increasing with the mass of the resonance and therefore the 
probability of the small size configuration being oc l/ry is decreasing with an increase of mass (for fixed 
S,L). Besides, the asymptotic freedom in QCD dictates decreasing of the coupling constant relevant for the 
behavior of the charmonium wave function at small relative distances. Experimentally one finds from the 
data on the leptonic decay widths that 1/fv drops by a large factor with increasing mass. An additional 
suppression arises due to the weakening of the soft exclusive nondiagonal VN ^ V'N amplitudes between 
states with the different number of nodes. Hence, it seems possible to determine the imaginary parts of 
diagonal rescattering amplitudes from analysis of the data using the GVDM equations. Since in the medium 
energy domain the energy dependence of soft rescattering amplitudes is well reproduced by a factor s'^"^, 
one can determine the real parts of the amplitudes using the Gribov-Migdal relation (Eq. [T^ . 

The above analyses demonstrate that the relative importance of the non-diagonal transitions increases 
an increase of the quark mass. It would be of interest to investigate experimentally the case of strangeness 
production. It is likely to correspond to e substantially larger than ep ~ 0.2 but much smaller than e,//^ ~ 1.7. 



3 Onset of perturbative color opacity at small x and onium co- 
herent photoproduction. 

Interaction of small size color singlet objects with hadrons is one of the most actively studied issues in 
high-energy QCD. The QCD factorization theorem for exclusive meson electroproduction at large Q^, and 
J/ip, T photoproduction 0310^1 allows to evaluate the amplitude of the production of a vector meson by a 
longitudinally polarized photon 7^ -I- T — > V + T through the convolution of the wave function of the meson 

^If the s-wavc mechanism dominates in the charmonium production in other hard processes the ratio would be a 

universal number. Since i/i" can be easily observed via its characteristic DD decays in a number of the lepton, hadron and 
nucleus induced processes the study of the discussed process can provide a new way to probing dynamics of charmonium 
production in various reactions. 
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at the zero transverse separation, hard interaction block and the generahzed (skewed) parton density The 
LT approximation differs strongly from the expectations based on Glauber model approaches and on the two 
gluon exchange models, because it accounts for the dominance of the electroproduction of spatially small 
quark-gluon wave package and its space-time evolution which leads to formation of a softer gluon field. The 
latter effect results effectively in an increase of the size of the dipole with increase of the energy. 

In perturbative QCD (similar to QED) the total cross section of the interaction of small systems with 
hadrons is proportional to the area occupied by color within projectile hadron 50 leading to the expectation 
of the color transparency phenomenon for various hard processes with nuclei. The cross sections of incoherent 
processes are expected to be proportional to the number of nucleons in the nuclei, while the coherent 
amplitude is proportional to number of nucleons times the nuclear form factor. Possibility to approximate 
projectile heavy quarkonium as colorless dipole of heavy quarks can be formally derived from QCD within 
the limit when mass of heavy quark mg — > oo but x — Attiq/s is fixed and not extremely small (36) . In this 
kinematics the size of heavy quarkonium is sufficiently small to justify applicability of PQCD. 

For practical purposes the crucial question is at what squeezing becomes effective. Probably the most 
sensitive indicator is the i-dcpendence of the meson production. The current HERA data are consistent 
with the prediction of 46, 36 that the slopes of the p and J/^p production amplitudes should converge to 
the same value. This indicates that at small x configurations much smaller than average configurations in 
light mesons {d ~ 0.6/m) dominate for p-meson production at > 5 GeV^, while the J/tp production is 
dominated by interaction in small size configurations for all Q^. Therefore, one expects the regime of color 
transparency for x > 0.03 where the gluon shadowing is very small/absent. 

Recently the color transparency (CT) phenomenon was observed at FNAL by E791 experiment |5T] 
which studied the coherent process of dissociation of a 500 GeV pion into two jets off the nuclei. The 
measurement has confirmed a number of predictions of |49) including the A-dependence, and the transverse 
and longitudinal momentum distributions of the jets. Previously the color transparency type behavior of 
the cross section was observed also in the coherent J/4> photoproduction at {E~^) = 120GeV [52j . 

A natural question is whether the color transparency will hold for arbitrary high energies? Two phenom- 
ena are expected to work against CT at high energies. One is the LT gluon shadowing. There are theoretical 
expectations (see discussion below) supported to some extent by the current analyses of the data on DIS 
scattering off nuclei (which do not extend deep enough into the shadowing region) that the gluon distri- 
butions are shadowed in nuclei as compared to the nucleon: Ga{x,Q'^)/AGn{x,Q^) < 1. This obviously 
should lead to a gradual but calculable in QCD disappearance of color transparency j^BJ 1311 , and to onset 
of a new regime, which we refer to as the color opacity regime (one can think of this regime also as a regime 
of generalized color transparency since a small qq dipole still couples to the gluon field of the target through 
a two gluon attachment and the amplitude is proportional to the generalized gluon density of the nucleus). 
Another mechanism for the violation of CT at high energies is the increase of the small dipole-nucleon cross 
section with energy cx Gn{x,Q'^). For sufficiently large energies this cross section becomes comparable to 
the meson-nucleon cross sections. One may expect that this would result in a significant suppression of 
the hard exclusive diffractive processes as compared to the LT approximation. However it seems that this 
phenomenon is beyond the kinematics achievable for the photoproduction of J/i/i-mesons in UPC of heavy 
ions at RHIC {x « 0.015, Q\fj « AGeV^) but could be important at heavy ion UPC at LHC. 

Hence, a systematic study of the onium production in the coherent scattering off nuclei at collider energies 
will be very interesting. One should emphasize here that with decrease of the size of the onium the eikonal 
(higher twist) contributions die out quickly (provided x is kept fixed). In particular, for the T case one 
probes nuclear gluon fields at the transverse scale of the order of 0.1 fm or Qlff ~ AOGeV^. The J/ip 
case is closer to the border line between the perturbative and nonperturbative domains. As a result the a 
nonperturbative region appears to give a significant contribution to the production amplitude |53| . 

Let us discuss the photoproduction amplitude j + A ^ J/ip{T) + A in more details. We are interested 
here in the Wjp range which can be probed at LHC. This region corresponds to rather small values of x. 
In this situation interaction of qq pair, which in the final state forms a quarkonium state is still rather far 
from the BEL. Hence the key problem in the theoretical treatment of the process is taking into account 
the nuclear shadowing. A number of mechanisms of coherent interactions with several were suggested for 

•^Proportionality of the hard diffractive amplitudes to the gluon density of the nucleon was discussed for hard pp diffraction 
in 1471 . and for J/ip production )48| in the BFKL approximation and for the pion diffraction into two jets )4H| in the leading 
log approximation 1491 . 
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this process. We focus here on the the leading twist mechanism of shadowing which is presented by 
diagrams in Fig^ There exists quahtative difference between the mechanism of interaction of a smaU dipole 
with several nucleons and the case of a similar interaction of an ordinary hadron. Let us for example consider 
interaction with two nucleons. The leading twist contribution is described by the diagrams where two gluons 
are attached to the dipole. To ensure that nucleus remains intact in such a process, the color singlet lines 
should be attached to both nucleons. These diagrams (especially the one of Fig. ^p) are closely related to 
the diagrams describing the gluon diffractive parton densities (which are measured at HERA) , and hence to 
the similar diagrams for the gluon nuclear shadowing 0]. 

The amplitude of high energy coherent heavy onium photoproduction is proportional to the generalized 
gluon density of the target, Gt{xi,X2, t, Q\j^), which depends on the light-cone fractions xi and X2 of gluons 
attached to the quark loop. They satisfy a relation: 

9 

Xi — X2 = = X (19) 

s 

If the quark Fermi motion and binding effects were negligible, X2 xi. The resolution scale Qe// — "^g' 
where nig is the mass of the heavy quark. Numerical estimates for the photoproduction of J/ip give Q^ff ~ 
3 — 4GeF^|2niinS| reflecting a relatively small mass of c-quark and indicating that this process is on the verge 
between nonperturbative and perturbative regimes. On the contrary, the mass of the beauty quark is huge 
on the scale of soft QCD. In this case hard physics dominates, attachments of more than two gluons to bb are 
negligible and the QCD factorization theorem provides a reliable tool for description of the T production. 
This is especially true for the ratio of the cross sections of T production from different targets since the 
higher twist effects due to the overlapping integral of the bb component of the photon and T cancelled out 
in the ratio. As a result, in the leading twist shadowing approximation the cross section of the process 
J A — > TA is proportional to the squared nuclear gluon density distribution and can be written in the form 

2 u„.„ 



(^-fA-*VA[S) = — 



GA{xi,X2,Qij:f,t = 0) 



AGNiXx,X2,Qlff,t ^ 0) 



dt 



d^dze'^'-^e-'^'pib^z) 



2 



(20) 



Numerical estimates using realistic potential model wave functions indicate that for J/ip, x\ ^ 1.5a;, X2 ~ 
a:/2 and for T, X2/X1 ~ 0.1 Modeling of the generalized parton distributions at moderate 

suggests that to a good approximation G{xi,X2,t = 0) can be approximated by the gluon density ai x — 
{xi -|-X2)/2 jJBJIM]. For large and small x GPDs are dominated by the evolution from Xi{init) ^ Xi. Since 
the evolution conserves xi — X2i effect of skewedness is determined primarily by the evolution from nearly 
diagonal distributions, see Ref. ^1] and references therein. In the case of the T production it increases 
the cross section by a factor ^ 2 [561 157] and, potentially, this could even obscure the connection of the 
discussed effect with the shadowing of the nuclear gluon densities. However, the analysis of shows that 
the ratio of GPD on a nucleus and on a nucleon at t = is a weak function of X2 , slowly dropping from its 
diagonal value (x2 = x\) with the decrease of X2- Overall this observation is in a agreement with the general 
trend mentioned above that it is more appropriate to do comparison of diagonal and non-diagonal cases at 

X = {xi+ X2)/2. 

In the case of double scattering contribution (Fig. QJa) there is another way to address the question of 
the accuracy of the substitution of the ratio of generalized gluon densities by the ratio of diagonal parton 
densities at the normalization scale. We notice that the amplitude corresponding to Fig. is expressed 
through the nondiagonal matrix element of the diffractive distribution function, {x]p,Q'^,xi,x,t), which is 
an analog of generalized parton distribution. In the diagonal limit of xi — X2 it coincides with the diffractive 
gluon distribution. It depends on the light-cone fraction which nucleon lost in \in) and {out\ states : xp 
xp — X respectively, f3in = xi/xp, Pout — (xi — x)/{xp — x), and t,Q^. If we make a natural assumption 
that 

g^{xi,x,xp,Ql,t) = \Jg^iPin,Ql,xp,t)g°{(3out,Ql,xp ~ x,t), (21) 

we find that numerically in the kinematics we discuss the resulting skewedness effects are small as compared 
to the uncertainties in the input gluon diagonal diffractive PDFs. 

Hence in the following we will approximate the ratio of generalized gluon densities in the nucleus and 
nucleon by the ratio of the gluon densities in nucleus and nucleon at a; = rriy/s. In the case of T use of 
i ~ a;/2 maybe more appropriate. This would lead to slightly larger shadowing effect. 
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It was demonstrated in 0] that one can express the quark and gluon nuclear shadowing for the interaction 
with two nucleons in a model independent way through the corresponding diffractive parton densities using 
the Gribov theory of inelastic shadowing [20] and the QCD factorization theorem for the hard diffraction 
|59| . An important discovery of HERA is that hard diffraction is indeed dominated by the leading twist 
contribution and gluons play a very important role in the diffraction(this is loosely referred to as gluon 
dominance of the Pomcron). Analysis of the HERA diffractive data indicates that in the gluon induced 
processes probability of the diffraction is much larger than in the quark induced processes 4 . The recent 
HI data on diffractive dijet production provide an additional confirmation of this observation. Large 
probability of diffraction in the gluon induced hard processes could be understood in the s-channcl language 
as formation of color octet dipolcs of rather large sizes which can diffractively scatter with a quite large cross 
section. The strength of this interaction can be quantified using optical theorem and introducing 



IGtt dadiff{x,Q'^,t„iin) _ 167r 



atot{x,Q^) dt (l+772)GAr(:r,g2) 



dx]pg^{ — :X]p,Q'^,tmin), (22) 

J Xp 



for the hard process of scattering of a virtual photon off the gluon field of the nucleon. Here rj is the ratio of 
the real to imaginary parts of the elementary diffractive amplitude, is the momentum scale determining 
virtuality of the gluons, xp is the momentum fraction of the pomeron with the corresponding cut-off scale 
Xp — 0.03, and g^i^^ i tmin) is the diffractive gluon density distribution of nucleon which is known 

from the HI analysis of the diffractive data at the scale « 4 GcV^. 

An important feature of this mechanism of coherent interaction is that it is practically absent for x > 
0.02 0.03 and may rather quickly become important with decrease of x. 

We calculated the ratio of the gluon density distributions in Eq. (201 using the leading twist shadowing 
model ^ (for the details of calculations see Ref. 61 ). As the first step the nuclear gluon density distribution 
with account of the leading twist shadowing is calculated at the starting evolution scale Qq = 4 GeV^ 
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(23) 



As an input we used the HI parameterization of ff^(-^, xp, Qq, t,nin)- The effective cross section aeff{x, Qq) 
accounts for the elastic rescattering of the produced diffractive state off the nuclear nucleon and is determined 
by Eq. |221 Numerically it is very large at the starting scale of the evolution, see Fig. [51 and corresponds to 
the probability of the gluon induced diffraction close to 50%. This indicates that at the initial scale of the 
evolution interactions in the gluon sector are close to the BBL for x < 10~^ for the nucleon case, and even 
more so for the nuclei, where similar regime should hold for a larger range of the impact parameters. 

Note that the double scattering term in Eq.[2Slfor the nuclear parton densities satisfies QCD evolution, but 
the higher order terms do not. That is if we use different starting scale of evolution we would obtain different 
results of G{x, Q^)- The reason is that the terms cx c"//: n>2 are sensitive to degree of fluctuations in the 
cross sections of interaction of diffracting states. These fluctuations increase with increase of . This effect 
is automatically included in the QCD evolution, and it leads to violation of the Glauber - like structure of the 
expression for the shadowing at > Qq. Approximation for the n > 3 rescattering of Eg .1231 corresponds to 
an assumption that fluctuations are small at Qq scale since this scale is close enough to the scale of the soft 
interactions, see discussion in 4 . Thus, at the second step of the calculation we use NLO QCD evolution 
equations to calculate the shadowing at larger Q^ using the calculation at Qq as a boundary condition. In 
this way we also take into account the contribution of the gluon enhancement at x 0.1 which influences 
the shadowing at larger Q^. Note here that proximity to the BBL in the gluon sector which is reflected in a 
large value of cTeff may result in corrections to the LT evolution which require further studies. 

First we calculate the ratios of the cross sections of coherent photoproduction of J/tp and T off nuclei and 
nucleon (Fig. [HJ . Such ratios do not depend on the uncertainties of the elementary cross sections and provide 
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Figure 5: The effective cross section sfiadowing in tire gluon channel, (Teff{x), at = 4 GeV^ as a function 
of the Bjorkcn x for HI parameterizations of the gluon diffractive density. 
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Figure 6: The a;-dependence of the ratio of J/i/' and T production off Ca and Pb in Glauber model to that 
in Impulse Approximation (I A). Calculation with the HI parameterization|60| of the diffractive PDF. 
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a sensitive test of the role of the LT shadowing effects. In a case of the J/ij] photoproduction the gluon 
virtuality scale is 3-4 GeV^ with a significant fraction of the amplitude due to smaller virtualities [201 IHSI ■ 
Hence we will take the gluon shadowing in the leading twist at = AGeV^. Taking a smaller value of 
would result in even larger shadowing effect. In calculations for the T we take Q^ff = 4:0GeV^, though the 
result is not very sensitive to precise value of Q^ff since the scaling violation for the gluon shadowing for 
these is rather small. We find that in spite of a small size of T, which essentially precludes higher twist 
shadowing effects up to very small x, the perturbative color opacity effect is quite appreciable. It is worth 
noting that the effective cross section of the rescattering in the eikonal model is determined by the cross 
section of dipole -nucleon interactions at the distances d ~ 0.25 — 0.3/m in the J/ip case and ~ 0.1/m in 
the T case. These cross sections are 10 — 15mb for J/tp and ^ 3mb for T case for x ~ 10~^ and a factor 
of 1.5-2 smaller for x ~ 10"'^, see Fig. 13 in Ref. They are much smaller than the cross sections which 
enter into calculation of the gluon shadowing in the LT mechanism. 

We also estimated the absolute cross section of onium photoproduction for a wide range of the photon 
energies. This is of interest for the planned measurements of the onium production in the ultraperipheral 
heavy ion collisions at LHC. The dependence of the momentum-integrated cross sections on the energy 
W^N = -s/s is presented in Fig. |7| . In the case of J/ip production the calculations are pretty straightforward 
as the accurate data are available from HERA. The situation is more complicated in the case of the photo- 
production of T. So far the information about the elementary j + N ^ T + N cross section is very limited. 
There is the only ZEUS and HI data on total cross section for the average energy of y/s « lOOGeT^. Thus, 
to calculate the forward photoproduction cross section we used a simple parameterization: 

da,^^vj,{s,t) ^ / s_Y'' . (24) 

at V So / 

where sq = 6400 GeV^, the slope parameter Br = 3.5 GeV^^^ is fixed basing on the analysis of the two 
gluon form factor in Ref. 62 , and the energy dependence follows from the calculations of Ref. of the 
photoproduction of T in the leading log approximation with an account for the skewedness of the partonic 
density distributions. This elementary cross section is normalized so that the total cross section is in 

4 LARGE MASS DIFFRACTION IN THE LEADING TWIST 
LIMIT 

It is well known that inelastic diffraction at small t gives information on the fluctuations of strength in 
the projectile-target interactions 29, 63 . Application of this logic to the hadron scattering off nuclei have 
allowed to explain the A-dependence of two measured diffractive channels in pA, ttA scattering assuming that 
it coincides with the A-dependence of the total cross section of the inelastic diffraction, and the absolute 
total cross section of pA diffraction (the data exist for two nuclei only) at the energies 200-400 GeV, see 
[SHI for the review and references. With an increase of energy, the total cross section of NN interaction 
increases and fluctuations in the elementary amplitude lead to much smaller fluctuations of the absorption 
in the scattering off a heavy nuclei. As a result, one can expect much weaker A-dependence of the diffractive 
cross section in particular, adiff{p + A — > X + A) cx A"-^^ at LHC energies P as compared to 

(Tdiffip + A ^ X + A) cx A" *" at fixed target energies. For large produced masses we can also understand 
this suppression using the t channel picture of the Pomeron exchanges as due to the stronger screening of 
the triple Pomeron exchange, for review and references see |65j. 

Diffraction in deep inelastic scattering corresponds to the transition of the (virtual) photon into its 
hadronic components leaving the nucleus intact. Hence it is similar to clastic hadron- nucleus scattering 
rather than inelastic diffractive hadron-nucleus scattering. 

It was demonstrated recently that in the UPC collisions at LHC it would be possible to study nuclear 
parton densities using hard charm and beauty production in 7 + A interactions |66j . Naturally one can also 
use these and similar processes to measure diffractive parton densities of nuclei. Since these quantities in 
the leading twist satisfy the factorization theorem we can analyze them on the basis of the analysis of the 
diffraction in DIS. 

There is a deep connection between shadowing and phenomenon of diffractive scattering off nuclei. The 
simplest way to investigate this connection is to apply the AGK cutting rules [HZ!- Several processes con- 
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Figure 7: The energy dependence of the coherent J /i}) and T photoproduction off Ca and Pb in the LT 
approximation. 
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tribute to diffraction on nuclei: (i) Coherent diffraction in wliich the nucleus remains intact, (ii) Break-up of 
the nucleus without production of hadrons in the nucleus fragmentation region, (iii) Rapidity gap events with 
hadron production in the nucleus fragmentation region. In Ref. we found that for x < 3 ■ 10"'^, > 4 
GeV^, the fraction of the DIS events with rapidity gaps reaches the value of about 30-40% for heavy nuclei, 
with a fraction of the events of type (iii) rapidly dropping with A. 

We can use the information on cTeff for quarks and gluons to estimate probability of diffraction for 
different hard triggers at the resolution scale ~ Q§. First we consider the dependence of the fraction of the 
events due to coherent diffraction and due to the break-up of the nucleus on the strength of the interaction, 
ciff, neglecting fluctuations of the interaction strength. We find that for the realistic values of cr^jy the 
probability of coherent diffraction is quite large but increases with a^f / very slowly and does not reach the 
asymptotic value of 1/2 even for very large values of a (the later feature reflects presence of a significant 
diffuse edge even in heavy nuclei), see Fig. |H1 Thus, it is not sensitive to the fluctuations of (Jeff- We also 
found that the ratio of diffraction with the nucleus break-up and with the nucleus remaining intact is small 
(10-20%) in a wide range of nuclei, and slowly increasing with increase of (Jeff, see Fig. |S1 Hence it would 
require high precision measurements to constrain the dynamics using <Tq-ei/o'tot ratios. 

Comparing the values of the fraction of the diffractive events for quark and gluon induced processes off 
heavy nuclei and proton, we find that the relative importance of the quark induced events is increasing. 
Therefore, the scaling violation at large /3 for the diffractive quark distribution in nuclei will be stronger for 
nuclei than for a proton. Another interesting effect is that for heavy nuclei only genuine elastic components 
can be produced (inelastic diffraction is zero). Hence, the soft contribution at Qq due to triple Pomeron 
exchange is strongly suppressed see e.g. As a result, nuclear diffractive parton distributions at small /3 
are strongly suppressed (by a factor cx A^^^) at Qg though this suppression will be less pronounced at large 

due to the QCD evolution. This will lead to breakdown of the universality of the (3 distributions as a 
function of A. 

Though the diffractive parton densities change rather slowly with leading to a weak variation of the 
diffractive cross sections with (modulus the scaling factor) the fraction of the diffractive events at fixed x 
should significantly drop with increase of due to a large increase of the inclusive nucleon parton densities 
and decrease of the nuclear shadowing. 

For example, let us consider ultraperipheral collisions (UPC) at LHC where one can measure a the process 
"f + A — > jeti +jet2 + X + A in the kinematics where direct photon process "f + g — + qq dominates. In this case 
if we consider the process at say pt ~ lOGeV/c corresponding to ^ lOOGeV^ the fraction of diffractive 
events will be of the order 10%. The background from the strong interaction originates from glancing 
collisions in which two nucleons interact via a double diffractive process pp pp + X where X contains jets. 
Probability of the hard processes with two gaps is very small at collider energies - even smaller than .the 
probability of the single diffractive hard processes, see e.g. [SHI- Therefore, we expect that the background 
conditions will be at least as good in the diffractive case as in the inclusive case considered in ^H]. Thus, 
it would be pretty straightforward to extract coherent diffraction by simply using anti-coincidence with the 
forward neutron detector, especially in the case of heavy nuclei, see discussion in jZO]- As a result it would be 
possible to measure in the UPC the nuclear diffractive parton distributions with a high statistical accuracy. 
It is important that in difference from the diffraction to a vector meson it would be possible to determine on 
the event by event basis the energy of the photon which induced the reaction, since the rapidity of the photon 
is close to the rapidity of two jets. As a result it would be possible to perform the measurements for large 
rapidities (selecting the events generated by a photon of higher of two energies allowed by the kinematics of 
production particles in the interval of rapidities yi < y < 1/2) and to determine diffractive parton densities 
for pretty small x. 



5 LARGE MASS DIFFRACTION IN THE BLACK BODY LIMIT 

One of the striking features of the BBL is the suppression of nondiagonal transitions in the photon interaction 
with heavy nuclei 17 . Indeed, in the BBL the dominant contribution to the coherent diffraction originates 
from "a shadow" of the fully absorptive interactions at impact parameters b < Ra, so the orthogonality 
argument is applicable. We use it to derive the BBL expression for the differential cross section of the 
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Figure 8: Dependence of the partial cross sections ratios for the hard processes on the effective cross section 
defined in Eq.2f compared with the shadowing of the inclusive hard process. 
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production of the invariant mass for scattering of (virtual) photons . For the real photon case: 

dtdM^ Stt 167r -tR\ ' ^ ' 

Here p{M^) = a{e^e~ hadrons) / a{e^ ji^ ^~). Comparison of the measured cross section of the 
diffractive production of states with certain masses with the BBL result (Eq. I25|l would allow to determine 
up to what masses in the photon wave function interaction remains black. Similar equation is valid in the 
BBL for the production of specific hadronic or quark-gluonic final states {qq, qqg, etc) in the case of the 
coherent nuclear recoil. This allows to measure any component of the light cone photon wave function 
which interacts with the BBL strength corresponding final states in the coherent processes. The onset of 
BBL limit for hard processes should reveal itself also in a faster increase with energy of cross sections of 
photoproduction of excited states as compared to the cross section for the ground state meson. It would 
be especially advantageous for these studies to use a set of nuclei - one in the medium range, like Ca, and 
another with A ^ 200. This would allow to remove the edge effects and use the length of about 10 fm of 
nuclear matter. 

One especially interesting channel is exclusive diffractive dijet production by real photons. One may 
expect that for the 'yA energies which will be available at EIC or at LHC in UPC the BBL in the scattering 
off heavy nuclei would be a good approximation for the masses M in the photon wave function up to few GeV. 
This is the domain which is described by perturbative QCD for x ~ 10~^ in the case of the proton targets, 
and larger x for scattering off nuclei. The condition of large longitudinal distances, a small longitudinal 
momentum transfer, will be applicable in this case up to quite large values of the produced diffractive mass. 
In the BBL the dominant channel of diffraction to large masses is production of two jets with the total cross 
section given by Ea. H25|l and with a characteristic angular distribution (1 + cos^ 6), where 6 is the cm. angle 
[2]. On the contrary, in the perturbative QCD limit the diffractive dijet production, except the charmed 
dijet production, is strongly suppressed |71l I72j . The suppression is due to the structure of coupling of the 
qq component of the real photon wave function to two gluons when calculated in the lowest order in as . As 
a result, in the real photon case hard diffraction involving light quarks is connected to production of qqg 
and higher states. Distribution of diffractively produced jets over invariant mass provides an important test 
of the onset of BBL limit. Really, in the DGLAP/CT regime differential cross section of forward diffractive 
dijet production should be oc and be dominated by charm jet production. This behavior is strikingly 

different from the BBL limit expressions of 7 . Thus, the dijet photoproduction should be very sensitive 
to the onset of the BBL regime. We want to draw attention that qq component of the photon light-cone 
wave function can be measured in three independent diffractive phenomena: in the BBL off the proton, in 
BBL off a heavy nuclei, in the CT regime where the wave function can be measured as a function of the 
interquark distance |49j . A competing process for photoproduction of dijets off heavy nuclei is production 
of dijets in 7 — 7 collisions where the second photon is provided by the Coulomb field of the nucleus. The 
dijets produced in this process have positive C-parity and hence this amplitude does not interfere with the 
amplitude of the dijet production in the jlP interaction which have negative C-parity. Our estimates indicate 
that this process will constitute a very small background over the wide range of energies . 

6 Coherent vector meson production in UPC at LHC 

Ultraperipheral collisions (UPC) of relativistic heavy ions at RHIC and LHC open a promising new avenue 
for experimental studies of the photon induced coherent and incoherent interactions with nuclei at high 
energies [ZHl ■ Really, the LHC heavy ion program 0] will allow studies of photon- proton and photon 
- nucleus collisions at the energies exceeding by far those available now at HERA for 7 — p scattering. 

Hence, we can analyze an opportunity to study the phenomena discussed above combining the theory of 
photo induced processes in the ultraperipheral AA collisions with our studies of the coherent photo(electro) 
production of vector mesons. We can use the standard Weizsacker- Williams approximation 75^ to calculate 
the cross section integrated over the momentum of the nucleus which emits the quasireal photons. 

The cross section of the vector meson production integrated over the transverse momenta of the nucleus 
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which emitted a photon can be written in the convoluted form: 



da{AA ^ VAA) 
dy 



Nj{y)<JjA^VA{y) + N^{-y)a^A^VA{-y)- 



(26) 



Here y is the rapidity 




(27) 



The flux of the equivalent photons (y) is given by a simple expression \7l^ : 



N{y)) 



^ I d'bTAA{b)^X'[K!{X) + iK',{X)] 



(28) 



Here Ko{X) and Ki{X) arc modified Bessel functions with argument X = ™^'^ , 7 is Lorentz factor and b 
is the impact parameter. The Glauber profile factor 



accounts for the inelastic strong interactions of the nuclei at impact parameters b < 2Ra and, hence, 
suppresses the corresponding contribution of the vector meson photoproduction. 

Recently the STAR collaboration released the first data on the cross section of the coherent p-meson 
production in gold-gold UPC at Wnn = V^nn = 130 GeV 76 . This provides a first opportunity to check 
the basic features of the theoretical models and main approximations which include the Weizsacker- Williams 
(WW) approximation for the spectrum of the equivalent photons, an approximate procedure for removing 
collisions at small impact parameters where nuclei interact strongly, and the model for the vector meson 
production in the jA interactions. In the case of the p-meson production the basic process is understood 
much better than for other photoproduction processes. Hence, checking the theory for this case is espe- 
cially important for proving that UPC could be used for learning new information about photon - nucleus 
interactions. Note here that the inelastic shadowing effects which start to contribute at high energies still 
remain a few percent correction at energies < 100 GeV relevant for the STAR kinematics. For LHC energy 
range one should account for the blackening of interaction with nuclei. In this case cross section of inelastic 
diffraction in hadron-nucleus collisions should tend to 0. So major impact for the calculation of the process 
of diffractive photoproduction of p meson would be necessity to neglect by the contribution of p' 

The calculated momentum transfer distributions at the rapidity y = and the momentum transfer 
integrated rapidity distribution for gold-gold UPC at y/s^N — 130 GeV are presented in Figs. Etjb jlOj . 
Let us briefly comment on our estimate of the incoherent p-meson production cross section. The momentum 
transfer distribution (dashed line in Fig.|51i) is practically flat in the discussed t± range. The total incoherent 
cross section obtained by integration over the wide range of t± is amc = 120 mb. To select the coherent 
production the cut t± < 0.02 GeV^ was used in the data analysis WEI. Correspondingly, the calculated 
incoherent cross section for this region of t± is amc — 14 mb. Our calculations of incoherent production 
which are based on accounting for only the single elementary diffractive collision obviously present the lower 
limit. The residual nucleus will be weakly excited and can evaporate only one-two neutrons. The events 
A + A ^ p + xn + Ai + A2 were detected by the STAR and identified as a two-stage process - coherent 
p-production with the subsequent electromagnetic excitation and neutron decay of the colliding nuclei |77|. 
In particular, the cross section estimated by the STAR for the case when only one of the nuclei is excited and 
emits several neutrons is a^^^ g„ = 95 ± 60 ± 25 mb. The momentum transfer distribution for these events is 
determined by the coherent production. Hence, it differs from that for incoherent events but in the region 
of very low t± it is hardly possible to separate them experimentally, and the measured cross section a!^^ g„ 
includes contribution of incoherent events on the level of 15%. 

The total rapidity-integrated cross section of coherent p-meson production calculated in the GVDM for 
the range of energies available at RHIC is shown in Fier. llOf solid line). We find cr*^^ = 540 mb at ^/snn — 130 
GeV. The value cr^^^ — 370 ± 170 ± 80 mb was obtained at this energy by the STAR from the data analysis 
at the low momentum transfer t± < 0.02 GeV'^. Thus, before making a comparison we should take into 
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Figure 9: (a) Momentum transfer dependence of the coherent and incoherent p-meson production in Au-Au 
UPC at y/sNN — 130 GeV calculated in GVDM. (b) Rapidity distributions for coherent p-meson production 
in the gold-gold UPC at ^^snn = 130 GeV calculated in GVDM. 
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Figure 10: Energy dependence of the total cross section for coherent p-meson production in the gold-gold 
UPC calculated in the GVDM. 
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Approximation 


Ca-Ca at LHC(7 = 3500) 


Pb-Pb at LHC(7 = 2700) 


Impulse 


0.6 mb 


70 mb 


Leading twist 


0.2 mb 


15 mb 



Table 1: Total cross sections of J/tp production in UPC at LHC. 



account this cut. It leads to a reduction of the cross section by « 10% (the dashed hue in Fig. ^j). In our 
calculations we didn't account for the t_L-dependence of the elementary amplitudes which are rather flat in 
the considered range of energies and momentum transfers as compared to that for the nucleus form factor. 
So, in the region of integration important for our analysis it is reasonable to neglect this slope. Nevertheless, 
an account of this effect would slightly reduce our estimate of the total cross section. Also we neglected 
a smearing due to the transverse momentum of photons and the interference of the production amplitudes 
from both nuclei '78' . 

This latter phenomenon results in the narrow dip in the coherent i_L-distribution a,t t± < 5 ■ 10^^ GeV^. 
All these effects do not influence noticeably the value of the t_L-integrated cross section but can be easy 
treated and taken into account in a more refined analysis. Thus we find cr'^^ = 490 mb to be compared to 
the STAR value a^^^ — 370 ± 170 ± 80 mb. Since our calculation does not have any free parameters, this 
can be considered as a reasonable agreement. 

It was suggested in ^ to look for color opacity phenomenon using J/ip (photo) electroproduction. 
This however requires energies much larger than those available at the fixed target facilities and would 
require use of electron-nucleus colliders. At the same time estimates of the counting rates performed within 
the framework of the FELIX study have demonstrated that the effective photon luminosities generated in 
peripheral heavy ion collisions at LHC would lead to significant rates of coherent photoproduction of vector 
mesons including T in reaction 

A + A + A + V. (30) 

As a result it would be possible to study at LHC photoproduction of vector mesons in Pb-Pb and Ca-Ca 
collisions at energies much higher than the range W^p < 17.3 GeV covered at the fixed target experiment at 
FNAL [S2]. Note that even current experiments at RHIC (W^p < 25 GeV) should also exceed this limit. As 
it is clearly indicated by the STAR study the coherent photoproduction, leaving the both interacting nuclei 
intact, can be reliably identified by using the veto triggering from the two-side Zero Degree Calorimeters 
which select the events not comprising the escaped neutrons. The additional requirement which enables to 
remove contribution of the incoherent events with the residual nucleus in the ground state is selection of 
the produced quarkonium with small transverse momentum. In Fig. II II we compare the momentum transfer 
distributions for the coherent J/ip and T photoproduction calculated in the Leading Twist shadowing model 
with the corresponding distributions for incoherent photoproduction. Note that we estimated the upper 
limit of incoherent cross section simply as the free elementary cross section on the nucleon target multiplied 
by the number of nucleons A. 

In Fig. I12b .b we present the rapidity distributions of the J/^ coherent production for peripheral collisions 
at LHC calculated including effects of gluon shadowing and in the impulse approximation. At the central 
rapidities we find suppression by a factor 4 for a case of Ca and more strong, by a factor 6 for Pb. The total 
cross sections are given in Table ^ The rapidity distributions for coherent T production in the UPC with 
Ca and Pb beams are shown in Fig. 1121;. d and the corresponding total cross sections are given in Table |21 
As it is seen from comparison of the Leading Twist shadowing based calculations to that performed in the 
Impulse Approximation the yield of T is expected to be suppressed by a factor 2 at central rapidities due to 
the leading twist nuclear shadowing. 

Hence, study of the coherent photoproduction of the heavy quarkonium states at LHC opens an important 



Approximation 


Ca-Ca at LHC (7 = 3500) 


Pb-Pb at LHC (7 = 2700) 


Impulse 


1.8 fib 


133 fib 


Leading twist 


1.2 fib 


78 fib 



Table 2: Total cross sections of T production in UPC at LHC. 
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Figure 11: The momentum transfer distribution for the coherent J/V' and T production in Ca-Ca and Pb-Pb 
in UPC at LHC. Solid Hne - calculation with the Leading twist shadowing, dashed line - the momentum 
transfer distribution for the incoherent photoproduction. 
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Figure 12: The rapidity distributions for the J/tp and T coherent production off Ca and Pb in UPC at LHC 
calculated with the leading twist shadowing based on HI parameterization of gluon density(solid line) and 
in the Impulse Approximation(dashed line) 
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avenue for investigating the nuclear gluon distributions and the shadowing effects in the kinematics which 
would be very hard to probe in any other experiments. 



7 CONCLUSIONS 

We demonstrated that coherent diffraction off nuclei provides an effective method of probing a possible onset 
of BBL regime in hard processes at small x. We predict a significant increase of the ratio of the yields of 
p, p' mesons in coherent processes off heavy nuclei due to the blackening of the soft QCD interactions in 
which fluctuations of the interaction strength are present. An account of nondiagonal transitions leads to 
a prediction of a significant enhancement of production of heavier diffractive states especially production of 
high pt dijets. Study of these channels may allow to get an important information on the onset of the black 
body limit in the diffraction of real photons. We argued that the fluctuations of strengths of interactions has 
been observed at intermediate energies in the diffractive photoproduction of vector mesons. We discuss the 
opportunity to look for the transition from the nuclear color transparency to the regime of the color opacity 
in the ultrarelativistic peripheral ion collisions at LHC. 

We thank J.Bjorken, A.Mueller, G.Shaw, members of the UPC study group, and our coauthors M.McDermott, 
A.Freund, V.Guzey, L.Gerland, for useful discussions and GIF and DOE for support. M.S. and M.Z. thank 
INT for hospitality while this work was completed. 
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